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We calculate the thermal evolution of tt — tt scattering lengths, in the frame of the Nambu-Jona- 
Lasinio model. The thermal corrections were calculated at the one loop level using Thermofield 
Dynamics. We present also results for the pion thermal mass. Our procedure implies the modeling 
of a propagating scalar meson as a resumation of chains of quark bubbles, which is presented 
explicitly. We compare our results with previous analysis of this problem in the frame of different 
theoretical approaches. 
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I. INTRODUCTION 

The Nambu-Jona-Lasinio model (NJL) is an effective, 
non renormalizable, low energy description of hadron dy- 
namics ([l|, 0]). The main motivation for this model is 
the natural way in which the phenomenon of chiral sym- 
metry breaking appears. The model has the same global 
symmetries of QCD, being therefore appropriate for de- 
scribing low energy hadronic phenomenology. Of special 
interest are processes involving pions, since they play a 
quite special role among hadrons. They are the lightest 
hadrons, being also pseudoscalar Goldstone bosons. 

In heavy ion relativistic collisions, it turns out that the 
pion multiplicity dominates by far any other hadronic 
signature, being dominant in the particle content of the 
central rapidity region, where we expect the quark gluon 
plasma to be created. This suggests that pion interac- 
tions, i.e. scattering processes, will be important. Since 
those scattering processes occur in the presence of a ther- 
mal environment, our discussion should be presented in 
the frame of finite temperature quantum field theory. 
The main parameters associated to low energy tt-tt scat- 
tering are the so-called scattering lengths, which depend 
on the isospin channel. The first measurements of tt- 
tt scattering lengths were done by Rosselet et al A 
review about the present experimental status of tt-tt scat- 
tering lengths can be found in 

In this paper we address once again the problem of 
determining the thermal behavior of the tt-tt scattering 
lengths, in the frame of the NJL model, using the Ther- 
mofield Dynamics (TFD) formalism. For details about 
this formalism see [ll , @ • The novelty of our discussion 
is mainly technical, including the resumation in TFD of 
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bubble diagrams associated to a propagating sigma me- 
son, how to regularize products of distributions, as well 
as other non trivial intermediate steps that appear in the 
calculation of the diagrams. 

It should be noticed that the thermal behavior of the 
TT-TT scattering lengths in the NJL model has being dis- 
cussed long time ago in the imaginary time formalism, or 
Matsubara approach 0]. Our results, specially for 02 (T), 
do not agree with this paper. We do agree qualitatively, 
however, with the prediction for the thermal behaviour 
of TT-TT scattering lengths in the frame of Chiral Pertur- 
bation Theory [8|, and with a recent analysis based on 
the Linear Sigma Model . Our results indicate that the 
variation of the scattering lengths as function of tempera- 
ture, in particular for 02 (T), begins at much lower values 
of temperature than the behavior found in and . 

The plan of this article is the following: after present- 
ing the general formalism, we classify the relevant dia- 
grams for getting the tt-tt scattering lengths according to 
the NJL model. Then, we proceed to extend our calcula- 
tion to the finite temperature scenario. For this purpose, 
a resumation procedure of quark bubbles is presented in 
order to describe the thermal sigma meson. As a byprod- 
uct of our analysis we found also the thermal evolution 
of the pion mass. Finally, we compare our results with 
previous articles about this subject. 

II. SCATTERING LENGTHS AND THE 
RELEVANT DIAGRAMS AT ZERO 
TEMPERATURE 

In general, the tt-tt scattering amplitude can be 
parametrized according to 

Ta!3;Si = ^(S, t, u)5a0Ss^ -f A{t, S, u)5a'y5ps 

A{u,t,s)5as5i3-i. (1) 
where the a, (3, 7, 8 denote isospin components. 



By using appropriate projection operators it is possible 
to get the following isospin dependent scattering ampli- 
tudes 

T" = 3A{s, t, u) + A{t, s, u) + A{u, t, s) (2) 
= A{t,s,u) - A{u,t,s) (3) 
T"^ = A{t,s,u)+ A{u,t,s), (4) 

where denotes a scattering amplitude in a given 
isospin channel. 

Below the inelastic threshold the partial scattering am- 
plitudes can be parametrized as [lO| . 

where 5i is a phase-shift in the I channel. In fact our last 
expression can be expanded according to 




The parameters a\ and b\ are the scattering lengths 
and scattering slopes, respectively. In general the scat- 
tering lengths satisfy |ao| > |ai| > \a2\ > ■■■■ If we are 
only interested in the scattering lengths Aq , it is enough to 
calculate the scattering amplitude in the static limit 
or at threshold i.e. when s 4to^, t and u 

a'o^^T' {.s^4ml,t^0,u^0). (7) 

Our Lagrangian is given by 

Cnjl = ? (i^ - mo) V + G ((^V)' + (?i75T^)') , (8) 

where -0 denotes the quark fields, G is a coupling constant 
and niQ denotes the current quark mass. We will also 
consider an effective coupling g^^qq between the internal 
quark lines and the external pions. Notice that the model 
does not include a field associated to the scalar sigma me- 
son. The propagator of the sigma meson is represented 
in the random phase approximation (RPA) trough a geo- 
metrical sum of quark bubbles's chains. See fig. [1] where 
we show the effective exchange of a sigma meson in the 
s— channel. 

The lowest order diagrams that contribute to the tt-tt 
scattering lengths in the NJL model are shown in fig. [21 
where the lines with an arrow inside the loops denote the 
quarks, and the double lines the sigma meson. The exter- 
nal legs are of course the physical pions with momentum 
p which couple to the quarks trough the effective coupling 
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FIG. 1: Effective meson exchange modeled as a chain of quark 
bubbles. 
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FIG. 2: Relevant Feynman diagrams for the calculation of the 
TV ~ TV scattering lengths. 



At zero temperature these diagrams have been calcu- 
lated by Schulze [ll|. 

It is important to remark that the fermion propaga- 
tors depend on the constituent quark mass to, instead 
of the current quark mass toq, because our calculations 
arc done in the chiral broken phase. The quarks acquire 
their constituent mass due to a condensation of quark- 
antiquark pairs that appear in a completely analogous 
way to the Cooper pairs in the BCS theory of super- 
conductivity. This quark-antiquark pairing mechanism 
implies the appearance of the constituent mass trough 
the so called gap equation 

1 = + 8tGN,Nf I , (9) 

where TV = NcNf, being Nc = 3 the numbers of colors 
and Nf the numbers of flavors. In our case we will take 

III. THERMOFIELD DYNAMICS AND THE 
EFFECTIVE SIGMA MESON PROPAGATOR. 

In TFD we have to double the number of degrees of 
freedom in order to express any thermal average in the 
ensemble as a vacuum expectation value in a populated 
vacuum. The extra fields are called thermal ghosts and, 
therefore, the propagators in general will be given by a 
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two by two matrix. In the case of fcrmions. the matrix 
propagator has the form 



Sf{1) 



Sii Si2 

S21 S22 



(10) 



where 



Sii 



(/ 



TO 



2tt5 {f 



le 



e|iol/T + i 





e\io\/T^l - 

i 2ttS {P - 



le 




(11) 



where e is the sign function. Note that the Sn element 
is exactly the fermion Dolan- Jackiw propagator [l2| . We 
may define Sii{p) — Sp{p) + Sp{p), separating exphcitly 
the zero temperature from the finite temperature contri- 
bution. The other components of the matrix propagator 
appear only in internal lines which do not couple to ex- 
ternal physical particles. Therefore, in our case, we will 
have to deal with the full matrix propagator structure, at 
the one loop level, only when diagrams associated to the 
exchange of a sigma meson are considered. As we said, 
it is given as a sum of chains of quark bubbles. Notice 
also that in general, there are two different coupling con- 
stants for the four quark vertices, depending on the type 
of quark lines that come together 01, 0- 

We will now construct the effective thermal sigma me- 
son propagator. In general, the chain of bubbles will 
imply a resumation of the type X]n=o where 



D 



Dn D12 
D21 D: 



22 



(12) 



is a matrix whose elements Dij represent a bubble with 
an "open vertex" of type i at the left side of the bubble, 
and with a vertex of type j at the right side = 1, 2). 
In order to avoid double counting of vertices when powers 
of D are considered, we have to introduce quark bubbles 
with an open vertex only at the left side. It is clear that 
each element D^j involves the sum of the two possibilities 
to build a chain with two bubbles, in such a way that we 
have an open vertex of type i at the left side and a vertex 
of type j at the right side. This means that the two pos- 
sible values for the vertex in between are included. The 
same happens for every matrix element of some power n 
of D. 

It is convenient to define the following matrix 



Z = 



1 



1 - D 
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£'22- 

Vl{l-D22) 
V2D21 



D11D22 
V1D12 

V2{l-Dn) 



D12D21 



(13) 



that essentially corresponds to X^^^o where, we have 
explicitly introduced the vertices Vi associated to the left 



side of the sum of the quark bubbles. A similar treatment 
was introduced previously in the discussion of thermal 
renormalons in the At/)'* theory [Tsj . 

In order to appreciate the role of the Z matrix let us 
consider diagrams (d), (e) and (f) of fig. [2l There we 
have quark triangles that couple to the effective sigma 
meson trough open vertices. To be more explicit, let us 
denote by Ti the triangle diagram with an internal open 
vertex of type i. Then we introduce the vectors 



= ( Ti r2 ) 



Y = 



ri 
r2 



(14) 



In this way, the complete amplitude for these kind of 
diagrams, including both the zero as well as the finite 
temperature parts, is given in this matrix language by 



T = Y^ZY. 



(15) 



In the case of the triangle diagrams (rf), (e) and (/) of 
fig. m we need the following elements 



in" + 2in^ I «i 

i i 



Dn = 

D12 = 
D21 = 



D22 = (in" - 2in^ ) V2 



(16) 

(17) 
(18) 

(19) 



to build the Z matrix. In the above expression 
-m°(-m'3) denotes a zero (finite) temperature single 
bubble contribution. Notice that D12 = D21 = be- 
cause there is no support for the corresponding integrals, 
as can be cheeked, since we are using constituent quark 
masses m, being to > 2to^. In order to handle other sin- 
gular contributions, like products of delta functions with 
the same argument, it was necessary to take different 
external pion momenta, i.e. a kinematical configuration 
away from the threshold. In this way, singularities dis- 
appear and we are able to compute the corresponding 
integrals, finding that they vanish identically when the 
limit to the threshold configuration is taken. Because of 
the same reasons, the triangle diagram T2 which couples 
to the sigma meson trough a vertex V2, vanishes. In this 
way, for the triangle diagrams we found 



r2 



4r 



/3 



where 



2to 



' d ^2 \i>-^qqJ 



A ;2 
dl- 



"~mE^-miy 



(20) 
(21) 



(22) 



where Ei = vP + n? and np is the usual Fermi-Dirac 
distribution 



npiz) 



1 



.\zo\/T 



(23) 
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FIG. 3: The single quark bubble that model a pion in the 
NJL approach. 



Wc have introduced a temperature dependent coupling 
gf^^q between pions and quarks. This coupling can be 
obtained as 



an 



ps 



Qp2 



(24) 



— (m^ ) 



where —iWpg is the single quark bubble (pseudoscalar 
case) shown in fig. [31 At this point it is worthwhile to 
remember that such a quark bubble is used to model a 
pion in the NJL approach. T selects the isospin channel: 
Ti = Tj ~ T3 corresponds to the 7r° and/or Ti = r^*^ 
and Tj = t'^^^ are associated to the 'k^'^^ . In the previous 
equation, is the thermal pion mass which is obtained 
in the next section. 



IV. TEMPERATURE CORRECTIONS TO THE 
TT-TT SCATTERING LENGTHS 

Following the notation introduced by Schulze [ll|, all 
diagrams at zero temperature are given in terms of the 
three following integrals 



K{p) 
L{p) 



dH 


1 


{2TTf {12 _ ,„2 _ 


Vie) [{l+pf -m^ + iej 


dH 


1 


(27r)^ (/2 _ ^2 _ 


- ie)^ {^l + p)^ — + iej 


dH 


1 


(27r)^ (;2 _ ^2 _ 


- ie)^ (^l + p)^ — + ie^ 



(25) 



where A is a 3-momentum cutoff. After performing the 
angular and the Iq integrals, we found 



- 2^ 



- 4^ 



dl 



Ei{AEf~p^) 
^ ^, e {p^ l2Ef) 
Ef{4Ef-p^f 
P {20Ef - p2) 



dl 



Ef{AEf-p^Y 



(26) 
(27) 
(28) 



In the NJL approach, the gap equation plays a fun- 
damental role. At finite temperature it acquires the 



1 = 



^ + 8^GN 



m 



dH 



HGN 



{2^f P - 
^ dH inpiEi) 



- le 



{2^r 



El 



(29) 



The thermal corrections to the effective sigma meson 
propagator, at the level of a single quark bubble are 
shown in the fig. [31 but where the 175 T has been replaced 
by the identity matrix according to the NJL Lagrangian. 
This loop is given by 



-Uf^ (p) = 2iN 
i 



d^l npiEi) 
(27r)^ El 



(4to 



where Jiji) is 




1 



dl 



PupjEi) 
El (p2 - 4£;2) 



(31) 



Notice that at this level the Dolan-Jackiw propagators 
are enough. In particular, the thermal contributions arise 
from the insertion of only one pure thermal propagator 
Sp in the loop. The other propagator has to be the 
normal one at zero temperature. It is easy to see that two 
thermal propagators do not contribute. The denominator 
Dz{p) that appears in the Z matrix is given by 



Dz{p) = l-2GIil^{p)^AGIil^{p). 



(32) 



Using the thermal mass gap equation (|29p this expression 
can be written as 

Dz{p) - — + AiGNp" {I (p) + 2iJ (p)) . (33) 

TO 

Following the same procedure, we could also construct 
an effective pion propagator, which is not necessary for 
our purpose of getting the tt-tt scattering lengths. From 
this analysis, however, we get the thermal pion mass de- 
termined from the pole of the effective propagator 



Too 



1 



m 4iGN {I (to^) + 2iJ (to^)) 



(34) 



The behavior of is shown in fig. [H This allows us to 
express the quotient toq/to in terms of the other quan- 
tities that appear in the previous equations. Replacing 
Toq/to in (|33)) we get 



Dzip) = AiGN (p^ - 4to2) (/ (p) + 2iJ (p)) 

-AiGNml {I {m„) + 2iJ (m^)) . (35) 

An expression for to^ was also found trough the mass gap 
equation in the Matsubara or imaginary time formalism 
by previous authors We agree qualitatively with 

them. 

It is important to remark that in equation p5p . m^r is 
the usual pion mass at zero temperature because we are 
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FIG. 4: Thermal evolution of the pion mass. 



implementing a calculation at the one loop level. The 
inclusion of in the last expression goes beyond this 
approximation. 

Finally, we get the following expression, consistent at 
the one loop perturbation theory, for the effective sigma 
meson propagator 



1 AiGN 
2ND zip)' 



(36) 



In addition to the triangle diagram given in eq. (|22|) . 
we need, also associated to the diagram [2{e), the follow- 
ing quantities 



^ ./O 



A ;2 



PnF{Ei){AE^ +ml) 
d fmPnF{Ei)Ei 



dm^ V (4£^f " ml) 



dl 



dl 



(37) 



where 



^ . PnF{Ei){mW + (4m^ + 2ml)Ef ~ 8Ef) 



Ei{4Ef-ml)^ 



For the diagrams (d), (e) and (/), fig. [21 we get 



(38) 



■0+/3 



i?^'' (2m.) (rS + 4r^) 

T^' = D"+^{0){T^^+T'l+2T'iy 



0+/3 



(39) 



where [lli 



-87Vm (7(0) - mlK{m^)) {g° 



qqi 



(40) 



(7°^^ can be read from the first term in eq. [24l In terms 
of our integrals I{p) and K{p), it can be rewritten as 



9nqq 



-Ni {l{0) + I{m^)~mlK{m^))) 



-1/2 



(41) 



We find g°^g ^ 3.56. 

Concerning the temperature corrections to the box di- 
agrams (a), (b) and (c) in fig. [2l we only need the Dolan- 
Jackiw fermion propagators. For the same reasons we 
mentioned previously, in each box diagram will survive 
only the insertion of one pure thermal fermion propaga- 
tor S^, which can be any of the four propagators inside 
the box. Two, three or four thermal insertions vanish 
identically. The results are 



n-O+P 
7-O+/3 



a 



(42) 



where 



Jo 



\2^i9^qqJ Q^. 

.4 



PupjEQEi 
(AEf-ml)^ 

_d_ f Pnp{Ei)Ei \ 
, (9m2 1^ (4£;f -ml)) 



^(AE?-mlV ■ 



7-/3 



^2 Ei{AE\-.,.^, 

/-^ d ( Pnp{Ei)Ei{AEt 
" i^Ef-ml)^ 

UEf) 



^2 \y-Tqq. 



dim? 

4 PnpiEi) (m: 



^-^\dl 



5?- 



m2)3 



(43) 



7^*^ and were also obtained in [Tl|, and are given by 
T," = 47V(.9^,,)'(m2i^(TO,)-/(0)-/(m.)) 

= 8iV {gPj' {2mlK{m^) - /(O) - \mlL{m,)^ . 

(44) 

Using these results, we are able to find the thermal tt — 
TT scattering lengths. For each diagram in figUl the total 
amplitude is given as the sum of the zero temperature 
and the finite temperature parts. In this way we get 



7-0+/3 

rrO+P 
7-O+/3 



0+/3 



7-0+/3 

e 



7-0+/3 

r," + 47;'3^ 

Dl+f' (2m.)(r3 + 4r 

Dl+i' (0) fr2 + rf 
7-0+/3 



2rf 



(45) 
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FIG. 5: Thermal dependence of the tt — vr scattering lengths, 
normalized to its zero temperature value. The upper curve 
correspond to a^'^ /a^, and the lower one to a^^^ /a2. Tem- 
peratures are in MeV. 

where the superscri pt refers to zero temperature terms, 
which are given in To express them, we have used 
our functions /, K and L. Finally, the scattering lengths 
are given by 

= ^{^%'^' + ^r^^')- (46) 
Their thermal dependence is shown in fig. [5l 



The parameters we used in the calculation are 
A = 631 MeV, m = 339 MeV, (= (T = 0)) ~ 
138 MeV. This set of parameters implies also mo — 5.5 
MeV (the average between the lightest quark masses, 
TO„ ~ 4 MeV, ma- 7 MeV), U = 93 MeV and G ~ 5.51 
GeV^^. These values for the parameters are standard in 
the literature for the case with a 3-dimensional cutoff 0] . 

At zero temperature we found Oq — 0.161, 02 ~ 
—0.043. These numbers are in agreement with Wein- 
berg's results HI = Imlli2-Kfl = 0.16 and af = 
— 2mJ/327r/^ = —0.044. However, they disagree with 
the experimental results [3], ao = 0.26 ± 0.05 and 02 = 
—0.028 ± 0.012. Nevertheless, our goal here was to find 
the thermal evolution of the scattering lengths normal- 
ized by the zero temperature values. 

Our results agree qualitatively well with two previous 
calculations. The first one Q, was done in the frame of 
Chiral Perturbation Theory at the one loop level, whereas 
the Linear Sigma Model was used in the second analy- 
sis 0. In our case, the scattering lengths vary faster as 
function of temperature compared with the other two pa- 
pers. However, we disagree qualitatively with the results 
obtained in [31 . This calculation was done also in the 
frame of the NJL model, but using the imaginary time 
formalism to compute the thermal corrections. 
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